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ABSTRACT 

Stars more massive than about 8 end their lives as a Supernova (SN), an event of fundamental importance Universe-wide. Theoret- 
ically, these stars have been expected to be either at the red supergiant, blue supergiant, or Wolf-Rayet stage before the explosion. We 
performed coupled stellar evolution and atmospheric modeling of stars with initial masses between 20 Mq and 120 Mq. We found that 
the 20 Mq and 25 Mq rotating models, before exploding as SN, have spectra that do not resemble any of the aforementioned classes 
of massive stars. Rather, they have remarkable similarities with rare, unstable massive stars known as Luminous Blue Variables 
(LBV). While observations show that some SNe seem to have had LBVs as progenitors, no theoretical model had yet predicted that 
a star could explode at this stage. Our models provide theoretical support for relatively low-luminosity LBVs exploding as SN in the 
framework of single stellar evolution. This is a significant shift in paradigm, meaning that a fraction of LBVs could be the end stage 
of massive star evolution, rather than a transitory evolutionary phase. We suggest that type lib SN could have LBV as progenitors, 
and a prime example could be SN 2008ax. 
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1. Introduction 

Stars more massive than about eight times the mass of the 
Sun, although rare, have been key players in the cosmic his- 
tory through their enormous input of ionizing photons, energy, 
momentum, and chemical species into the Universe. Their short 
lives end with a core-collapse Supernova (SN), a generally lumi- 
nous event that can be traced up to cosmological distances. 

The morphological appearance of massive stars before the 
SN event is very uncertain. Observationally, red supergiants 
(RSG) have been confirmed as SN progenitors for stars with 
up to 18 Mq, while Wolf-Rayet (WR) stars have not yet been 
confirmed as SN progenitors (Smartt 2009). Theoretically, pro- 
genitors of core-collapse SN have been typically related either to 
the c lass of WR. blue supergiant (BSG). or RSG stars ( Georgv 



ttie c lass 01 WK. blue supi 
et aL l2012tlLangedl2012l) . 



Luminous Blue Variables (LBV), a 
heterogeneous class of unstable massive stars, have also been 
observationally connected to SN progenitors. This link has been 
made base d on the infer ence of wind variability from the radio 
-lightcurve dKotak & yinb.2006) and line profiles (iTrundle et al.l 
llOOl iGroh & Vinkll2011h. presence of a dense circumstellar 
medium (Smith et al."2007'), and photometry of the progeni- 
tor (Gal- Yam & Leonard 2009; Mauerhan et al. 2012). How- 
ever, this has not been supported by stellar evolution theory yet, 
which had instead suggested LBVs to be in a transitory phase 
from being an O-type star burning hydrogen in its core to becom- 
ing a He-core burning WR star (Ekstrom et al. 2012; Maeder & 
Meynet 2000). Such a glaring discrepancy between observations 
and theory exposes a main gap in the understanding of massive 
stars, their deaths, and consequently their impact in the Universe. 

So far, stellar evolution models have been able to self- 
consistently predict the stellar parameters such as temperature 
(r*) and luminosity (L*) only up to the stellar hydrostatic sur- 



face. However, especially at the pre-SN stage, this layer is not 
reached directly by the observations because of the presence of 
an atmosphere and dense stellar wind. To properly classify the 
SN progenitor from stellar evolution models, we present here 
a novel approach, coupling radiative transfer atmospheric/wind 
calculations done with CMFGEN (Hil lier & Millerl il998) with 
thos e performed with the Geneva stellar evolution code ( Ekstrom 
et al.|2012|). Using as input the stellar structure conditions, we 
are able to generate a synthetic spectrum that can be directly con- 
fronted to the observations. As we discuss in this Letter, having 
a spectrum allows us to gain enormous insight into the nature 
of the progenitor Here we report a quite surprising theoretical 
finding that single stars with initial mass in the range 20-25 Mq 
have spectra similar to LBVs before exploding as SNe. 



2. Physics of the models 

We performed coupled stellar evolution and atmospheric mod- 
eling of stars with initial masses between 20 Mq and 120 Mq. 
Hereafter we focus on the results of the 20 Mq and 25 Mq ro- 
tating models. The analysis of the full range of masses will be 

reported elsewhere (Groh et al. 2013, in prepara tion). 

Th e evolutionary models are those from Ekstrom et al.l 
(1201 2h . computed with the Geneva stellar evolution code. They 
assume solar metallicity and initial rotational speed of 40% of 
the critical velocity for break-up. The prescription for the ro- 
tational diffusion coefficients is taken from IZahnI (11992) and 
iMaeded (IT997.) . The radiative mass loss rates are implemented 
following the pr escription frorn JVink et all (1200 lb . and for the 
RSG phase from lde Jager et al] (11988 ). Because of variations in 
the ionization level of hydrogen beneath the surface of the star 
during the RSG phase, some layers might exceed the Edding- 
ton limit, possibly driving instabilities. In this case, the radiative 
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mass loss is increased by a fac tor of three in our mod els, which 
matches the observations from lvan Loon et al. I (l200l . 

To compute the output spectra we used the spherically sym- 
metric, full line blanketed, non-local thermodynamical equilib- 
rium, atmosph eric radiative transfer code CMFG EN (Hillier & 
Miller 1998). For the hydrodynamics of the wind, we assume a 
standard yS-type law with /3 - 1, while a hydrostatic solution is 
computed for the subsonic portion of the atmosphere. The two 
solutions merge at 0.75 of the sonic speed. The wind terminal 
velocity is computed using the scaling with the escape velocity 
from Kudritzki & Puis (2000). Wind clumping is included with 
a volume filling factor (/) approach, and all models computed 
here assume / = 0.1. 

We use the outputs from the Geneva stellar structure calcula- 
tions, such as the radius, luminosity, mass, and surface abun- 
dances, as input in CMFGEN. For consistency, we adopt the 
same mass-los s rate recipe as t hat used by the Geneva evolu- 
tion code (from l^nk et al.ll2001l in the parameter range analyzed 
here). We use the temperature structure of the stellar envelope to 
merge the CMFGEN and the Geneva stellar structure solutions . 
Changing the location of the merging point between the CMF- 
GEN and Geneva solutions would change the temperature at the 
hydrostatic surface, which would make the SN progenitor look 
hotter or colder However, the conclusions reached here are in- 
dependent of this approach. For reasonable changes in the merg- 
ing layer, no significant difference in spectral morphology is ob- 
tained, with the pre-SN model spectra still resembling LB Vs. 



3. The fate of rotating 20-25 Mq stars: LBVs 
exploding as Type II SN 

Figure [T] shows the evolutionary tracks of rotating massive stars 
with initial masses of 20 Mq and 25 Mq at solar metallicity. At 
the pre-SN stage, the 20 Mq model has = 20 358 K and = 
1.9 X 10^ Lo while the 25 Mq has = 24625 K and = 
3.2 X 10^ Lq (blue filled squares in FiglU. Several classes of 
massive stars, such as BSG and LBVs, are found in this range of 
Ti, and L*. In the absence of a spectrum, conventional wisdom 
would have suggested a BSG (based on and L*) or a late-type 

(WN classification ( WNL, using the chemical abundance criteria 
from lGeorgvl2012h . 
Strikingly, we found that the pre-SN synthetic spectra of our 
models are remarkably similar to the observed spectra of LBVs, 
such as AG Carinae, HR C ari nae, P Cygni , HDE 316285, amon g 
others (iGroh et al.ll2009allbt iHillier et al. 1 119981: lNaiarroll200lh . 
Our results show that the spectral morphology of the 20 Mq and 
25 Mq pre-SN models is markedly different than that of RSG, 
BSG, or WR stars. Figure |2] shows the pre-SN synthetic optical 
spectra computed with our atmospheric modeling. To illustrate 
the morphology of a typical observed LBV spectrum, we also 
display the prototypical LBV AG Carinae in Fig. |2] 

The similarities are conspicuous, especially because our 

(models do not have initial mass, rotation rate, mass-loss rate, 
or wind terminal velocity tuned to reproduce the spectrum of 
AG Car or any other LBV. One can clearly see that the main 
features present in the optical spectrum of LBVs, such as stron 
H, He I, and Nn lines with P-Cygni profiles (iGroh et al.ll2009 
IStahl et anil993h . are also present in the synthetic spectra at the 
pre-SN stage. While some LBVs may display a WNl Ih spectral 
type at visual minimum, we note that He ii /14686 is not present 
in emission and, therefore, a WR classification with the latest 
pos sible spectral types (W N9h - WNllh) would be inappropri- 
ate (ICrowther et al.lll995h . We obtained that the morphological 
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Fig. 1. Evolutionary tracks of rotating models with initial masses 
of 20 Mo and 25 Mq at solar m etallicity (Z=0 .014), shown in the HR 
diagram and calculated in Ekstr om et al The initial rotation rate 

is 40% of the critical velocity for break-up. Blue filled squares indicate 
the pre-SN stage, and representative evolutionary phases are indicated. 
Temperatures refer to those at the last envelope layer. 



appearance of the SN progenitor is regulated not only by and 
Li, but also by the stellar wind properties, such as the mass-loss 
rate and wind terminal velocity. In particular, the modeled stars 
have relatively high mass-loss rate (M ~ 1 -4x 10"''Moyr"') and 
low wind terminal velocity (voc ~ 270 to 320 km s"'), which are 
in line with values typically encountered in LBVs. The model 
properties at the pre-SN stage are summarized in Table [T] 

The Galactic LBVs span a range of luminosities, with 
log(L*/L0) ^ 5.2 - 6.5 (Clark et al. 2005). Observationally, it is 
unclear whether low- and high-luminosity LBVs represent well- 
defined, separated classes. Our results seem to indicate that low- 
luminosity LBVs {\og{Li, I Lq) ^ 5.2-5.5) are the endpoints of 
stellar evolution, while high-luminosity LBVs could still evolve 
to the WR phase. 

In addition to spectroscopic similarities, both the 20 Mq and 
25 Mq rotating models are extremely close to the Eddington 
limit, a limit at which the luminosity of the star is so high that 
the gravity could not longer hold the star in hydrostatic equilib- 
rium. Our modeling indicates that the Eddington parameter V 
is around 0.85 - 0.90 in the deep subsonic layers of the wind 
(at a Rosseland optical depth of 20). This occurs because of the 
relatively high luminosity-to-mass ratio resulting from the evo- 
lution, which in turn supports the prescribed mass-loss rates and 
wind-terminal velocities used in the stellar evolution and spec- 
tral calculations. We recall that the proximity to (or surpassing 
of) the Eddington limit seems to be a key factor in the LBV phe- 
nomenon (Grafener et al. 2012; Groh et al. 2011; Humphreys 
& Davidson 1994; Smith & Owockill2006h . Our assumed mass- 
loss rates may be affected by the distinct properties of stellar 
mass loss near t he Eddington limit dGrafener & Hamannll200l 
IVink et al.l201 11 . which could change the extension of the atmo- 
sphere and the effective temperature. However, for reasonable 
changes, qualitatively we expect similar results, i. e., a spectral 
type reminiscent of LBVs would still be seen. 

A third kinship between our models and LBVs are the chem- 
ical abundances at the surface. Both models show He and N 
overabundances and H, C, and O depletion, which is a signature 
that CNO-processed mateiial is present on the suiface (Fig. |3]l. 
For the 20 Mq model, the He and N abundances (0.74 and 0.0075 
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Fig. 2. Synthetic optical spectra of the 20 Mq (black) and 25 Mg (blue) rotating models at the pre-SN stage. The flux has been continuum- 
normalized and offset for clarity. The strongest spectral features are identified. We also show the ob served s pectrum of the prototypical LBV AG 
Carinae obtained in 2002 July 04 (green), when the star had an effective temperature of 16400 K dOroh e t al. 2009b). Note that AG Car has a 
higher luminosity than the pre-SN models, and is shown only to illustrate the morphology of a typical LBV spectrum. We note the remarkable 
similarity between our synthetic spectra and a typical LBV spectrum. This suggests that some LBV stars are caught just before the explosion. 



by mass fra ction, respectively ) are remarkably similar to those of 
LBVs (e. g. lClarket an2009l) . The 25 Mq model predicts a more 
chemically enriched surface with He and N (0.92 and 0.016 by 
mass fraction, respectively) and, as expected, its spectrum shows 
stronger He I and N II lines than that of the 20 Mq model or the 
LBV observations. Still, the 25 Mq model spectrum looks more 
similar to LBVs than to any other class of massive stars. 

We note that the LBV classification is phenomenological, 
where one of the following properties has to be present. Either 
the star shows S-Doradus type variability, or giant eruptions of 
sev eral solar masses must have occurred ( Humphrevs & David- 
sonl l994tlvan Genderenll2001h . Therefore, had the SN model 
progenitors been observed, they would technically be classified 
as "candidate LBVs", since none of the above phenomena are in- 
cluded in our calculations. Despite semantics, here we show that 
the model SN progenitors are alike stars in the quite rare LBV 
phase in many aspects, such as spectral appearance, proximity to 
the Eddington limit, and chemical abundances at the surface. 

One may wonder why the rotating model explodes as an 
LBV, while the coiTesponding non-rotating model 20 Mq does 
not evolve back to the blue and explodes as a RSG, showing a 
type IIP SN. This is mainly an effect of rotational mixing which 
produces a larger He-core. The He-core is the inner part of the 
star, enclosed in the first layer encountered going deeper into the 
star where the He mass fraction is greater than ~ 0.75. It is a 
well known effect that the larger is the mass fraction of the total 
mass occupied by the He-core, the bluer is the position of the 
star in the HR diagram (Giannone 1967). In our rotating 20 Mq 
model, the He-core encompasses 100% of the total mass, while 
in the non-rotating corresponding model, it encompasses only 
71%. Hence the difference of position in the HR diagram. 



4. Implications for massive star evolution 

Our results indicate that for rotating massive stars with initial 
masses between 20-25 Mq, a previously unrepoited evolution- 
ary sequence occurs. The star appears at the Zero-Age Main Se- 
quence as late-type O star (06V to 07.5V), evolving to a BSG/O 
supergiant (OSG), a RSG, a BSG/blue hypergiant (BHG, for 
the 20 Mq model) or O supergiant (OSG)/WNL (for the 25 Mq 
model), and finally appearing as an LBV before the SN event. 



20 Mq: 07.5V BSC RSG -» BSG/BHG -> LBV SN; 

25 Mq: 06V -> OSG RSG -» OSG/WNL LBV SN. 
Our models indicate that the pre-mortem LBV phase is short 
(~ 5000 years ), which is consistent with the o bserved rarity 
of LBVs (Clar k et al.ll2005L iHumphrevs & Davidson 1994; van 
Genderen 200 ij" The 20 Mq and 25 Mq models lose most of 
their mass at the RSG phase (8.7 Mq and 9.6 Mq respectively), 
which enables the star to have the aforementioned surface prop- 
erties required to show an LBV spectrum at pre-death. Thus, 
mass loss at t he RSG phase has a key impact in the subsequent 
evolution (see lGeorgvl20 1 2h . 

Let us look at the expected properties of the SN arising from 
the explosion of these low-luminosity LBVs. For this purpose, 
we investigate the chemical composition of the star before the 
explosion. Figure [3] shows the chemical structure of our 20 Mq 
rotating model at the end of the core C burning phase. The star 
consists of a He-rich core containing 94% of the total mass and 
extending over slightly more than 10% of the total radius of the 
star, suiTounded by a very shallow H-rich envelope containing 
about 6% of the mass and extending over about 90% of the total 
radius. We see that very little amount of H remains at the surface. 

According to lGeorgv et al.l (l20I2h . the 20 Mq model presents 
a very nice fit with the observed charac teristics of the core- 
collapse SN 2008ax (ICrockett et al.ll2008l) . These are (a) a CO 
core mass between 4 and 5 Mq (the model gives 4.73 Mq), (b) a 
small amount of hydrogen in the ejecta (the model gives 0.02 Mq 
of H) and (c) the estimated posi tion of the progenit or in the color- 
magnitude diagram (Fig. 10 of lGeorgv etanl20T2 ). Our revised 
colors for the rotating 20 Mq model (B - V = -0.14) confirms 
the very good agreement. SN 2008ax was found to be a type lib, 
i.e. an event which transitions from type II to lb. This suggests 
that type lib SN could have an LBV as progenitor, and a prime 
example could be SN 2008ax. 

Thus, for the rotating 20 Mq model, one would expect that 
the SN progenitor, which has an LBV-type spectrum, will ex- 
plode as a type lib SN. The 25 Mq model presents a very similar 
structure as the one shown in Fig.|3] thus will likely explode as a 
type lib SN too. Because the progenitor star has significant mass 
loss during the RSG phas e and crosses t he infamous "Yellow 
Void" in the HR diagram (Ide JageiilI99"8h . enhanced mass loss 
may have occuiTed in both phases. As a consequence, this would 
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Table 1. Properties of the SN progenitors from rotating models at solar metalicity, which have an LBV spectrum and explode as SN Il-b / Il-n. 
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He (surface) 


(Mo) 


(Mo) 


(Myr) 


(10* Lo) 


(K) 


(K) 


(Rg) 


(Re) 


(Moyr-') 


(kms-1) 


(mag) 


(mag) 


(mass fraction) 


20.0 


7.1 


10.5 


1.9 


20355 


19540 


35.3 


38.1 


1.2 X 10^^ 


272 


-6.8 


-0.14 


0.74 


25.0 


9.6 


8.6 


3.2 


25790 


20000 


28.4 


47.0 


4.6 X 10-5 


326 


-7.2 


-0.11 


0.92 



Notes. From lEkstrom et all ( 120 12h . Defined here as the temperature (radius) of the last layer of the stellar envelope. Defined here as the 
temperature (radius) of the layer where the Rosseland optical depth is equal to 2/3. 



produce a H-rich circumstellar medium with dense, slow veloc- 
ity material, and its interaction with the SN blast wave could 
produce a type Iln SN instead of a type lib. 

The (still scarce) observational evidence that suggests that 
LBVs can be SN progenitors indicates that the progenitors could 
come from a range of initial masses. Some of the LBVs sug- 
gested as SN progenitors h ad their initial mas s poorly con- 
strained, such as SN 2005gj (iTrundle et al.ll2008l) . while others 
were Unked to very massive stars (> 50 Mq), such as SN 2005g l 
^Gal-Yam & Leonard'2009^ and SN 2006gy (ISmith et al.ll2"007[) . 
Interestingly, some SNe have been proposed to have LBVs of 
modest initial masses as progenitors, compara ble to the mo dels 
we compute here, as in the case of SN 1987A (ISmithll2007t) . 

While our models still do not support LBVs as progeni- 
tors of SNe from stars with initial mass above 25 Mq, they 
provide unexpected theoretical support for explosions of low- 
luminosity LBVs in a single star scenario. The number of low- 
luminosity LBVs that explode as SN depends on the exact mass 
limits and distribution of rotational velocities at birth, since the 
non-rotating models explode as a RSG (20 Mq) or WNL/LBV 
(25 Mq). We note that the time average rotational velocities of 
the 20 and 25 Mq models during the Main Sequence (MS) are 
217 and 209 km s"'. This is comparable to the usual 200 km 
s ' considered as the mean velocities of such stars during the 
MS phase. Thus, one would expect that a significant fraction of 
stars in this mass range would follow the evolution previously 
discussed and explode as an LBV. 

In summary, our key theoretical finding is that stars with 
spectrum similar to LBVs immediately precede core collapse 
and the subsequent SN event of rotating stars with initial mass 
between 20 Mq and 25 Mq. Our models predict that the low- 
luminosity LB V at thepre-SN stage are blue, with B-V = -0.14 
to -0. 1 1 . They also have absolute magnitudes in the V-band be- 
tween My - -6.8 and -7.2. These values are three magnitudes 
fainter than those of SN progenitors that have been proposed to 
be LBVs, as in SN 2005h1 a nd 2009ip (My ~ -10: Gal- Yam & 
Leonard 2009; Mauerhan et al.l l2012h . Given the obvious obser- 
vational bias towards detecting the more luminous progenitors, 
a population of SNe with low-luminosity LBVs as progenitors 
could be awaiting detection, and the previous evidence for LBVs 
being SN progenitors could have revealed just the tip of the ice- 
berg. 
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